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ABSTRACT. 

An eigenspace assignment approach to the design of parameter insensitive control laws for 
linear multivariable systems is presented. The control design scheme utilizes flexibility in 
eigenvector assignments to reduce control system sensitivity to changes in system parameters. The 
methods involve use of the singular value decomposition to provide an exact description of allowable 
eigenvectors in terms of a minimum number of design parameters. In a design example, the 
methods are applied to the problem of symmetric flutter suppression in an aeroelastic vehicle. In this 
example the flutter mode is sensitive to changes in dynamic pressure and eigenspace methods are 
used to enhance the performance of a stabilizing minimum energy/ linear quadratic regulator 
controller and associated observer. Results indicate that the methods provide feedback control laws 
that make stability of the nominal closed loop systems insensitive to changes in dynamic pressure. 


tPartially supported by the grant ONR-N00014-86-K-0693 and by NASA Contract 
No.NASl-18107 while the author was in residence at ICASE, NASA Langley Research Center, 
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AN APPLICATION OF EIGENSPACE METHODS 
TO SYMMETRIC FLUTTER SUPPRESSION 


1. INTRODUCTION. 

A common problem in flight control applications is to design feedback and observer gain 
matrices so that the system remains stable and satisifies given performance requirements over a range 
of flight conditions. A basic approach is to fix a flight condition as a design point and to synthesize 
gain matrices so that stability and performance requirements are met at the design point.Nearby 
flight conditions yield evalution models. If performance requirements are not met at the evaluation 
points, one returns to the design model and adjusts the existing design to meet more stringent 
requirements. In such cases it frequently occurs that certain variables used for performance 
evaluation are more sensitive to changes in the design point than others and a reduction in sensitivity 
can be explicitly considered in the design process. In this paper eigenspace methods are presented 
for the design of parameter insensitive control laws. The methods are applied to the problem of 
symmetric flutter suppression in an aeroelastic vehicle. 

The ability to shape the fundamental modes of a system by modifying selected eigenvalues and 
eigenvectors underlies the appeal of eigenspace methods for control system design. The basic 
freedoms and limitations of eigenstructure assigment methods by full state feedback or output 
feedback are presented in the papers by Moore [12], by Srinathkumar [14], by Andry, Shapiro and 
Chung [2], and by Kautsky, Nichols and Van Dooren [7], The utility of eigenstructure assignment 
methods in observer designs is discussed in the work of Kazerooni and Houpt [8]. In these methods 
feedback and observer gain matrices are determined to yield a desired eigenstructure. Stability and 
transmission considerations motivate the location of desired eigenvalues, whereas the ability to shape 
the system response [2, 12], to enhance system performance [ 5], to reduce system sensitivity [15], 
or to design robust control laws [7, 8] motivates the selection of desired eigenvectors. 

In this paper, vector space methods are presented for eigenspace assignment which take into 
consideration eigenvalue sensitivity to plant parameter variations and performance constraints. The 
procedure is formulated as a constrained optimization problem in Section 2 along with a description 
of eigenspace assignment methods for the determination of feedback and observer gain matrices. 
Formulation of the procedure as a constrained optimization problem allows explicit consideration of 
costs associated with control effort in the design procedure. Eigenvalue sensitivity can be expressed 
explicitly in terms of closed loop eigenvectors and, consequently, the freedom available in assigning 
eigenvectors can be directly related to sensitivity reduction. In Section 2, the singular value 
decomposition [9] is used to provide a basis for the attainable eigenvectors associated with a desired 
eigenvalue. Thus an explicit parameterization/coordinatization of the attainable eigenvectors is 
obtained. An alternate algorithm for computing a basis for the attainable eigenvectors was presented 
in the work of Porter and DAzzo [ 13] and, more recently, by Kautsky, Nichols and Van Dooren 
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[7]. Some eigenvector assignment procedures require the control designer to specify a desired 
eigenvector and the nearest attainable eigenvector is then computed [2], [5]. The approach 
advocated, here, involves a precise description of the attainable eigenvectors and results in the 
display of the design freedom in terms of a minimum number of independent parameters. 

A flutter suppresssion/gust load alleviation problem is formulated in Section 3 and the 
results of an extensive design example are presented in Section 4. In this example, the flutter mode 
is sensitive to changes in dynmaic pressure and eigenspace methods are used to enhance the 
performance properties of a LQR/LTR designed compensator. Results indicate that the methods 
provide feedback control laws which make the stability of the nominal closed loop system 
insensitive to changes in dynamic pressure. 


2. PARAMETER INSENSITIVE CONTROL SYSTEM DESIGN. 

In this section a procedure based upon eigenstructure assignment methods for the design of 
parameter insensitive control laws is presented. A review of eigenstructure assignment methods 
useful in the design of full state feedback control laws and state observers is included. In this paper 
the singular value decomposition is used to provide a coordinatization of the allowable eigenvectors, 
which arise in full state feedback designs. These coordinates become the design parameters for 
sensitivity reduction. 

Eigenstructure assignment . Consider the linear system 

dx/dt = Ax + Bu+rrj (1) 


where A, B, and T denote matrices of appropriate dimensions and x, u, tj denote state, control, and 
disturbance variables respectively. Assuming full state feedback control, i. e. u = -Kx , the 
eigenvalues of A - BK can be arbitrarily assigned through a proper choice of the gain matrix K if 
and only if the system (1) is completely controllable. If (1) is uncontrollable then uncontrollable 
eigenvalues cannot be altered by a state variable feedback control law u = -Kx. 

The ability to assign eigenvalues and eigenvectors through state variable feedback control laws 
is summarized as follows [3, 7, 12, 14]. Let A = { A.^ : i=l, ..., n } be a self conjugate set of 
distinct complex numbers. Here A denotes a set of desired closed loop eigenvalues and must include 
all uncontrollable eigenvalues of the system. There exists a feedback gain matrix K such that (A - 
BK) Vj = A-i Vj for i = 1 , .... n if and only if 


i) ( v j }j_j n are linearly independent in C" andv. = v. if A.. = L 


ii) 


v. 


w. 


belongs to the Ker [ A^ I - A, B] for i = 1, ..., n 


= Kv. 

1 


for i = 1, ..., 


iii) 


w. 


n 
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In this case the feedback gain matrix K must satisfy K = [ wj, .... w n ] [ vj, v n ]'*. 
Thus A, is an eigenvalue for A - BK with corresponding eigenvector v provided 


[XI - A, B] 



(2) 


If W denotes a matrix whose columns form a basis for the Ker [ X I - A, B ] then 


v 

Kv 


= Wc 



c 


(3) 


where c denotes a vector whose components are the coordinates of [ v', (Kv)’ ]' with respect to the 
columns of W. Thus (3) provides a parameterization or coordinatization of the allowable 
eigenvectors. If A, is a controllable eigenvalue it may be shown that the dimension of 
Ker [ X I - A, B ] equals the number of columns of B and the columns of Wj , see (3) above, are 
independent if B is of full rank [12]. The singular value decomposition [7, 9], SVD, provides an 
efficient numerical method to compute an orthonormal basis for Ker [ X I - A, B ] . If the SVD of [ 
X I - A, B ] is 


[XI - A, B ] 


[U,,u 2 ] 


z 

0 


0 

0 




,v 0 ] 


(4) 


then the columns of W = V 2 form an orthonormal basis for the Ker [ X I - A, B ] . 

Equation (3) provides a parameterization or coordinatization of the attainable eigenvectors 
associated with X and, consequently, completely describes the freedom in assigning a particular 
eigenvector. One should note that the number of free parameters describing an atttainable 
eigenvector equals the rank( Wj) < dim Ker [ X I - A, B ] . These free parameters will be utilized 
as control design parameters in the application to follow. Further characterizations of the attainable 
eigenvectors and the number of degrees of freedom in eigenvector assignment are presented in the 
work of Kautsky, Nichols and Van Dooren [7]. 

In this design example the coordinates of an attainable eigenvector are modified to obtain 
feedback control laws which are insensitive to changes in a model parameter. The approach requires 
that derivatives of eigenvalues with respect to system parameters be computed. It is well known [7, 
10] that if A. is a distinct eigenvalue of a matrix Q(q), depending on a parameter q, and u’, v are left 
and right eigenvectors corresponding to A, , respectively, then 

dA. _ u’dQ/dqv 

dq u' v (5) 

The dependence of this sensitivity on the eigenvectors u’ and v indicates a potential reduction in 
sensitivity by a suitable choice of eigenvectors. Bounds on eigenvalue sensitivities are directly 
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related to the condition number of the modal matrix for the closed loop systems and are discussed 
in [7]. 

Ei gen structure assignment in observer desig ns. In a recent paper [8], Kazerooni and Houpt 
presented a procedure for loop transfer recovery based upon eigenstructure assignment of observers. 
This procedure for observer design will be used in the example to follow. These same methods have 
been previously applied in the design of an active flutter suppression system by Garrard, Liebst, and 
Farm [ 5 ]. For completeness the procedure is reviewed in this section. 

Suppose now only output feedback is available and an observer is to be employed to estimate 
the state, then equation (1) is coupled with the following output, observer, and feedback equations: 


y = Cx 

(6) 

dz/dt = Az + Bu + H(y-Cz) 

(7) 

u = -Kz 

(8) 


Here K and H denote feedback and observer gain matrices, respectively. 

Setting u = -K z + u 0 in place of (8) and e = x - z, one obtains 

de/dt = ( A - HC) e + B u 0 + rq . (9) 


In terms of transfer functions T zu (s) =T xu (s) - T eu (s). For recovery it is desireable to have 
T zu (s) = T xu (s) or T eu (s) = 0. 

The condition T eu (s) = 0 can be described in terms of transmission properties of (1) and (7). 
A complex number ^ is called an invariant zero [1 1] of the system (1) and (6) with left zero direction 
[V , p' ] provided 


[ V, p’ ] 


$1- A B 
C 0 


0 


( 10 ) 


and [v' , p’ ] * 0. If a gain matrix H can be chosen which assigns the eigenvalues {Xj : i=l, ...,n) 
of A - HC to invariant zeros { ^ j : i=l, ...» n) of (1) and (6) and left eigenvector to 
corresponding left zero direction Vj' then, necessarily. 


[VV.H] 


XA - A B 
C 0 


= 0 


Moreover, since the zero state response of (9) is 

V 1 V 

e(t) = 2^ u j v j’ Be * u o( l ) 

j=i 

where uj vj denote right and left eigenvectors of A corresponding to Xj, it follows that e(t) = 0, i.e. 
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T eu = 0. Thus for recovery it is desirable to place the eigenvalues of A - HC at the invariant zeros 
of (1) and (6) and the left eigenvectors at corresponding left zero directions. 

Numerical solution of (10) reduces to the solution of a generalized eigenvalue problem and 
existing computer software [5, 6, 7] can be used to obtain the invariant zeros and left zero directions. 
In this procedure the desired eigenstructure represents an ideal. There is no guarantee that the 
desired eigenstructure is attainable. In the design example to follow nearest attainable eigenvectors 
are computed using the methods presented in [6]. Once desired eigenvalues and attainable 
eigenvectors have been specified, the observer gain matrix can be computed using eigenstructure 
assignment methods. For further discussion on the application of this procedure see [ 3] and [6]. If 
the observer gain H(p) depends on a parameter p and H(p)/p -» BW as p — » °° where W is a 
nonsingular matrix then a standard arguement [10] implies that the eigenvalues of A - HC converge 
to the transmission zeros of (1) and (6). The control transfer function for the system described by 
equations (1), (6), (7), and (8) is u = - K(s)y with 

K(s) = K ( si - A + BK + HC ) _1 H. 

Let G(s) denote the transfer function associated with (1) and (7), i.e. 

G(s) = C ( si - A) _1 B. 

In this case Doyle and Stein [6] have also shown that K(s)G(s) — » KG(s) pointwise. Thus, from 
another viewpoint, it is desireable to place the eigenvalues of the observer at the transmission zeros 
of the system (1) and (6). 

Parameter insensitive control design . Here the objective is to determine the design parameters which 
minimize the sensitivity of certain performance variables while, at the same time, maintaining other 
performance variables within prescribed bounds. By introducing an appropriate penalty function, 
this problem can be formulated as a constrained optimization problem. 

In this paper, the design variables are those eigenvalues and eigenvectors of the system 
matricies A(q 0 ) - B( q c )K and A( q c ) - HC( q e ) which may be modified in order to reduce 
sensitivity. The design parameters are the real and imaginary parts of the designated eigenvalues 
and the coordinates of the associated eigenvectors. In the design process employed herein, a basis 
for the allowable eigenvectors must be computed for each eigenvector to be modified, see equation 
(3). The coordinates of an eigenvector with respect to this basis become the design parameters that 
relate to eigenvector selection. If a desired eigenvalue is to be left unaltered throughout a design, that 
is eigensystem design freedoms are not fully utilized, then this basis need only be computed once. 

Let a denote the vector of design parameters. In the design procedure certain performance 
variables, those which measure sensitivity to plant variations, are to be minimized. Let s( a) denote 
a vector of such variables. In the design example presented in Section 4, s( a) will denote 
eigenvalue sensitivity to plant parameter variation. Let p( a) denote a vector of performance 
variables to be kept within prescribed bounds. The prescribed lower and upper bounds for the i ^ 
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component of p( a) are denoted by LBj and UBj, i.e. it is desired that LBj < p( a)j < UBj. Let 
p(a) be the vector whose components are defined by p(a) = max{ 0, pfoOj-UBj, LB i -p(a) i }. 

The design procedure is to choose a to minimize the performance function 
J = 1/2 s( a ) ’ Q : s( a ) + 1/2 p( a ) Q 2 p( a ) 

where Qj and Q 2 denote diagonal weight matrices. Minimization of 1/2 s( a)' Qj s( a) tends to 

reduce sensitivity while the term 1/2 p (a)’Q 2 p (a) represents a penalty on performance constraint 

violations. The designer is free to vary the weighting matrices Q j and Q 2 in order to obtain an 
acceptable design. 

3. MODEL AND CONTROL PROBLEM DESCRIPTION. 

This work was motivated by the desire to attain flutter suppression and gust load alleviation in 
an aeroelastic vehicle. The planform of the vehicle wing with three control surfaces and three vertical 
accelerometers is depicted in Figure 3.1. State space equations for control design were obtained by 



Figure 3.1 Hypothetical Wing. 
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employing a rational s-plane approximation of the unsteady aerodynamic forces [1]. Data for the 
design and evaluation models was supplied by NASA. The state space equations are of the form 
dx/dt = Ax + Bu + rT] 
y = Cx 

where A is n x n , B is n x r, T is nxl and C is r x n with n = 26 and r = 3. The states x are 
associated with generalized positions (5), rates (5) and unsteady aerodynamic forces (5), control 
surface positions (3), rates (3), and accelerations (3), and gusts (2); the variables u and Tj denote 
control inputs and white noise into the Dryden turbulence model, respectively; and the observations 
are associated with vertical accelerometer measurements. The structure of the system matrices is 
depicted in Figure 3.2. The locus of system eigenvalues with variation in dynamic pressure, q, is 
depicted in Figure 3.3. Flutter onset occurs when the damping in the elastic mode becomes zero. 
The point chosen for control design is qq = 4.768 lb/in 2 . Other points serve as evaluation points. 
The control objective is to design feedback and observer gain matrices which stabilize the system 
over a range of values of dynamic pressure and which keep control surface deflections and rates 
within prescribed bounds. The maximum deflections and rates are set at 15° and 740° /sec, 
respectively. The control law is to be designed so that saturation of controls does not occur for an 
input gust spectrum having an rms intensity of 12 ft/sec. 
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Figure 3.2. Control System Structure. 

4. DESIGN EXAMPLE. 

At the design point qQ = 4.768 lb/in 2 , a stabilizing, minimum energy, full state feedback 
control law u = -K( qQ) x was determined using linear quadratic regulator methodology, i. e. 
weighting matrices for control and state variables were chosen to be the identity and zero matrices, 
respectively, see Table A1 of Appendix A for a listing of the gain matrix., With such a control law, 
unstable open loop poles are moved to their mirror image in the left half plane while other poles 
remain unaltered. 
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Figure 3.3 Locus of Open Loop Poles with Dynamic Pressure. 

At the design point the controlled system 
dx/dt = (A(q 0 )-B(q 0 )K(q 0 ))x + r(q 0 )Tl 

exhibits good performance and disturbance rejection properties. The rms values listed in column 
LQR of Table 4.1 indicate acceptable control surface activity and loads for a 12 ft./sec. rms gust 
velocity. In addition, the minimum singular value for the return difference matrix I + K( qj 
(sI-ACq^^^BCqo) over the frequency range 10 < co£ 300 is approximately one, indicating good 
disturbance rejection properties. 

At this design point an observer was determined by attempting to place the eigenvalues of A 
- HC at the transmission zeros of the plant, C( q 0 )( si - A( q c ))' 1 B( q c ), and corresponding 
eigenvectors at the left zero directions. The open loop eigenvalues, invariant zeros and actual 
assignments are listed in Table 4.2. Here eigenvalues 1-10 correspond to position and rate 
variables, 11-15 correspond to unsteady aerodynamic states, 16-24 correspond to actuator states, 
and 25-26 correspond to gust states. The zeros fall in three categories: those with large 
magnitude , those of the same order of magnitude and those with small magnitude relative to the 
magni tude, of the open loop eigenvalues. The six small magnitude zeros result from computattional 
errors and should be set to zero. Actual pole-zero assignments are listed in columns I and II of 
Table 4.2. The desired eigenvalues and eigenvectors for the assignment listed in column I of Table 
4.2 were made as follows: zeros of the same order of magnitude as those of the open loop poles 
(and corresponding left zero direction) were assigned to the nearest open loop eigenvalue (and 
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Control 

Deflections (deg) 


Control 

Rates (deg/sec) 


Bending Moment (in lbs) 
Shear (lbs) 

Torque (in lbs) 


LOR 

LQR/LTR 

.351 

.367 

.314 

.350 

.438 

.471 

41.9 

42.8 

36.9 

37.0 

51.6 

52.1 

26,344 

26,396 

492.0 

492.9 

1657. 

1708. 


Table 4.1 RMS Performance at Design Point for LQR and LQR/LTR Compensators. 



Eigenvalue 

Zero 

I 

II 

1 

9 

-13.8+430j 

-6.017+346.9j 

-6.017+346.9j 

-6.017+346.9j 

At 

3 

- 19.2+27 5j 

.1075 

-.12 

-4.8925 

4 

— 

-.1075 

-.1 

-5.1075 

5 

f. 

-1.4+194J 

-1.040+198.2j 

-1.040+198.2j 

-1.040+198.2j 

O 

7 

-43.6+136j 

.03479 

-.04 

-4.9652 

8 

— 

-.03479 

-.02 

-5.0348 

9 

10 

5.2+1 19.7j 

.000003+. 17j 

-.000003+. 17j 

-5.000003+.17j 

1U 

li 

-92 

-111.2 

-111.2 

-111.2 

12 

-113.6 

-129.8 

-129.8 

-129.8 

13 

-112 

-122.2+.06059j 

-122.2+.06059j 

-122.2+.06059j 

14 

-112.4 

-122.2-.06050J 

- 122.2- .06059j 

- 122.2- .06059J 

15 

-112.2 

-111.7 

-111.7 

-111.7 

16 

-179.9 

oo 

-879.9 

-879.9 

17 

-125.6+287.9j 

oo 

-825.6+287.9j 

-825.6+287.9j 

18 

— ■ 

oo 

— 

- 

19 

-175 

oo 

-875 

-875 

20 

- 127. 4+29 1.9j 

oo 

-827.4+29 1.9j 

-827.4+29 1.9j 

21 

— 

oo 

— 

- 

22 

-185 

oo 

-885 

-885 

23 

- 123. 9+283. 9j 

oo 

-823.9+283.9j 

-823.9+283.9j 

24 

— 

oo 

— 

- 

25 

! -.492 

-.492 

-.492 

-492 

26 

-.497 

-.497 

-.497 

-497 


Table 4.2 Eigenvalue-Zero Assignments for Observer Design. 
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eigenvector); zeros with small relative magnitude (and corresponding left zero direction) were 
assigned somewhat arbitrarily to the nearest unassigned eigenvalue (and eigenvector); the 
remaining eigenvalues were assigned a large negative real part and the corresponding eigenvectors 
were left in the open loop configuration. Note that in all cases the assigned eigenvalues have 
negative real part guarantying the stability of the observer. The resulting gain matrix H( q^) is 
listed in Table A3 of Appendix A. 

The system with observer is described by the equations 

dx/dt = A( q 0 ) x + B( q 0 ) u + T( q,,) Tl y = C( q G )x 

dz/dt = A( q^ z + B( q e )u +H( q^ ( y - C( q c )z ) u = - K( q c )z. 

The performance of the system with compensator is indicated by the rms values listed in Table 4.1, 
column LQR/LTR. Plots of the maximum and minimum singular values of the return difference 
matrix I + K(s)G(s) appear in Figure 4.1. 

Although the system with compensator appears to have good stability and performance 
properties the synthesized control laws are sensitive to changes in dynamic pressure. The locus of 
eigenvalue locations of A(q) - B(q) K(q 0 ) with variation in dynamic pressure is depicted in Figure 
4.2. Note that at the evaluation points q= 5.141, 5.537, 6.639 lb/in 2 the system matrices A(q) - 
B(q)K(q 0 ) are unstable. Additional calculations show that the matrices A(q) - H( q <) )C possess 


1.8 


1.4 


1 


.6 


.2 


Singular Values 


i 


10 100 Frequency 

Figure 4.1. Maximum and Minimum Singular Values of the Return 

Difference Matrix I + K(sl - A + BK + HC)' 1 HC(sI - A) _1 B. 
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Figure 4.2 Locus of Pole Locations for LQR Compensators A(q) - B(q)K(q ) . 

o 


unstable eigenvalues for q= 5.141 and 6.639 lb/in 2. It is desired to reduce this sensitivity while 
maintaining performance requirements by use of the eigenspace assignment methods discussed in 
the previous section. 

To illustrate the design procedure, a design involving only a small number of parameters is 
undertaken. The minimum energy stabilizing controller K( q 0 ) modifies the eigenstructrure of 
A( q Q ) in a simple manner, here the unstable eigenvalues are flipped symmetrically with respect to 
the imaginary axis while all other eigenvalue/eigenvector pairs remain unchanged. That is the 
stabilizing, minimum energy feedback control law only alters the mode shape associated with the 
open loop instability. In the design example, the closed loop eigenvalues are left in the stable 
location achieved by the minimum energy controller K( q 0 ) and an attempt is made to reduce the 
sensitivity of this control law to changes in dynamic pressure by modifying the eigenvector 
corresponding to the unstable mode. Letting X denote the desired eigenvalue associated with the 
unstable pole, the attainable eigenvectors v and modified gain matrix Kjyj must satisfy 


[ X I - A, B ] 


v 

k m v 


0 


In the example the dimension of Ker[X. I-A, B] is three. Since this is a complex vector space there 
are actually six free parameters which describe the attainable eigenvectors and modified gain matrix 
Thus the objective is to use the flexibility of eigenvector assignment to design a modified full state 
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feedback control law Kj^ which reduces sensitivity to changes in dynamic pressure, maintains 
performance requirements and is robust with respect to model uncertainties. Although all of these 
design objectives will not be met, the example does illustrate that many of them can be achieved by 
simply modifying the mode shape associated with the unstable eigenvalues. 

The sensitivity of the unstable eigenvalue to changes in dynamic pressure is calculated from 
the formula 

d X /dq = {u* (d( A(q) - B(q)K(q c )) /dq) v } / (u’ v) 

where u’ and v denote left and right eigenvectors of A(q) - B(q)K(q 0 ) corresponding to the 
unstable eigenvalue and d( A(q) - B(q)K(q 0 » /dq, is dX( q 0 ) /dq is estimated by the difference 
quotient (A(q) - B(q)K(q c ) - (A(q o ) - B( q Q )K(q 0 )}/(q - q^. The sensitivity is dX /dq = 43.88 + 
18.67j. The design procedure is to utilize the freedom in eigenvector assignment in order to 
minimize the magnitude of dX( qj /dq subject to the stability and performance constraints at this 
value of dynamic pressure. Robustness constraints are not explicitly taken into consideration. 

Let W = [ Wj’, W 2 ’ ]’ be a basis for the ker[ XI - A, B ] where X = -5.2 + 119.7j. The 
unstable open loop eigenvalue is 5.2 ± 1 19.7j. The attainable eigenvectors are linear combinations 
of the columns of Wj, i.e. v = Wj c for an arbitrary vector c. The design parameters a = c are 
the components of an attainable eigenvector with respect to the basis Wj, which is obtained using 
the SVD as described in Section 2. Here s( a) = d X/ dq. The components of the performance 
vector p( a) are a stability indicator, position and rate rms values for each controller, and wing 
bending moment, shear, and torque rms values. The stability indicator was taken taken as the 
maximum of the real part of the closed loop system eigenvalues. In our design example the 
prescribed lower bounds were LB = (-10,000,0,0,0,0,0,0,0,0,0) and the prescribed upper bounds 
were UB = (0,3.7,15,15,372,372,372,30000,1000,2000). The weight matrices were Q x = 1 and 
Q 2 = diag(10000,50,50,50, 25,50,50, 50, 50,50). 

A search over the six dimensional parameter space determines a desired eigenvector to 
achieve this minimu m and the corresponding gain matrix Kj^j. The sensitivity achieved at this point 
is dX( q Q ) /dq = 10.16 + 2.07j. The corresponding gain matrix is listed in Table A2 of Appendix A. 
To further illustrate the sensitivity reduction of this design the closed loop poles for the systems 
A(q) - B(q)Kj^j for various values of dynamic pressure have been depicted in Figure 4.3. Note that 
stability is achieved at each evaluation point Corresponding rms values are listed in Table 4.3. It is 
most interesting that these results can be achieved by simply modifying the mode shape 
corresponding to the unstable mode. The maximum and minimum singular values of the return 
difference matrix I + Kj^j C( q 0 ) (si - A( q c ) )' ! B( q G ) over a specified frequency range have 
been graphed in Figure 4.4. This figure indicates that the robustness properties of the minimum 
energy controller have not been maintained and that a more reasonable design would include 
singular value constraints. 
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Figure 4.3 Locus of Pole Locations for Modified LQR Compensators A(q) - B(q)Kjyf 



10 100 Frequency 

Figure 4.4. Maximum and Minimum Singular Values of the Return 
Difference Matrix I + K M (sI - A) _1 B. 
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Although the separation principle guarantees the stability of the compensator with feedback 
gain matrix K = and observer gain matrix H = H( q^ , instabilities still arise at off design 
points due to use of the observer gain matrix H = H( q^. Consequently a modified observer design 
was attempted, the desired pole locations for this modified design are listed in column II of Table 
4.2 and the resulting gain matrix is listed in Table A4 of Appendix A . The overall 
performance of the compensator may now be evaluated. First, it should be noted that numerical 
computations verify the stability of the compensator at all evaluation points. Thus the design 
procedure, i.e. modification of the eigenvectors associated with the unstable mode at the design 
point for feedback gain matrix design and assignment of observer eigenstructure to corresponding 
system zeros and invariant directions for observer gain matrix design has resulted in a 
compensator whose stability is less sensitive to changes in dynamic pressure. However one should 
also note that the rms values for the system with compensator, see Table 4.4, indicate a high level 
of control activity at the evaluation points furthest from the design point. Further improvement in 
performance should be attainable though modification of other mode shapes and by adjusting the 
weight matrices in the design procedure to obtain an acceptable blend between sensitivity reduction 
and control activity. Plots of the maximum and minimum singular values of the return difference 
matrix I + (si - A + BKjyj + Hj^C)'^Hj^C(sI - A)'*B are depicted in Figure 4.5. Again, 
robustness properties of the minimum energy controller have not been retained. 



Table 4.3 RMS Performance for Modified LQR Compensator dx/dt = (A(q)-B(q)Kj v j)x. 
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* saturation 

Table 4.4 RMS Performance for Modified LQR/LTR Compensator 

dx/dt = A(q)x-B(q)K M z+r(q)ri, dz/dt = H M C(q)x+(A(q)-B(q)K M -H M C(q))z. 



Figure 4.5. Maximum and Minimum Singular Values of the Return 

Difference Matrix I + K M (si - A + BK M + H M Q'^QsI-A^B . 
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CONCLUDING REMARKS. 

A procedure based upon eigenstructure assignment to reduce control system sensitivity to 
changes in system parameters has been introduced. Through an extensive design example it has 
been shown that modification of a single mode shape can lead to a gain matrix whose sensitivity to 
parameter variation is significandy reduced. The design example also indicates that a proper blend 
must be achieved between sensitivity reduction and performance constraints. Another problem that 
needs to be addressed is the explicit incorporation of robustness requirements into the design 
procedure. Current effort involves the incorporation of explicit robustness constraints in the 
design. Additional freedoms in eigenvector assignment are being used to achieve this objective. 
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Table A2. Modified Gain Matrix K, 
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Table A4. Transpose of Modified Observer Gain H, 
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